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A basic linearity of quantum dynamics, that density matrices are mapped hnearly 
to density matrices, is proved very simply for a system that does not interact with 
\Q . anything else. It is assumed that at each time the physical quantities and states 
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' are described by the usual linear structures of quantum mechanics. Beyond that, 
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■ the proof assumes only that the dynamics does not depend on anything outside 

the system but must allow the system to be described as part of a larger system. 
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that pure states are mapped one-to-one onto pure states and that entropy does not 
decrease. 
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I. INTRODUCTION 

Can we prove that quantum dynamics must be hnear? That the Schrodinger equation 
must be linear? What assumptions are needed? Is there a reason in principle that quantum 
dynamics is linear? If not, quantum dynamics could be a linear approximation of a nonlinear 
theory; then experiments might reveal small nonlinear effects. If linearity can be proved from 
first principles, most physicists will not consider nonlinearity as an option for theories or 
spend time and money on experiments to look for nonlinear effects. 

A definition of basic linearity for quantum dynamics is that density matrices are maped 
linearly to density matrices in an interval of time. This basic linearity is proved very simply 
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here for a system that does not interact with anything else. It is assumed that at each time 
the physical quantities and states are described by the usual linear structures of quantum 
mechanics. What is proved is the basic linearity for the changes of states and mean values in 
time. Beyond the structure at each time, the proof assumes only that the dynamics does not 
depend on anything that happens outside the system, or on the state of any other system, 
but must allow the system to be described as part of a larger system. The latter condition is 
close in spirit to the idea of envarianceQ], but the proof does not use quantum entanglement. 
The proof does not even require quantum mechanics; it could be done in classical mechanics 
as well. 

This proof of basic linearity is the first of two steps. It is done in Section |nl The second 
step, taken in Section IIIH links this basic linearity with previously established results to 
complete a simple derivation of the linear Schrodinger equation. For this it is assumed that 
the dynamics maps density matrices one-to-one onto density matrices in an interval of time. 
An alternative is to assume that pure states are mapped one-to-one onto pure states and 
that entropy does not decrease. 

The question of the linearity of quantum dynamics has a rich history. Wigner provedj^,!^ 
that quantum dynamics must be linear if it does not change absolute values of inner prod- 
ucts of state vectors. Bialynicki-Birula and Mycielskiljl proposed a nonlinear Schrodinger 
equation that inspired precise experimental tests^, 0,1^. Weinberg^, 0] proposeda general 
nonlinear form of quantum mechanicSjWhich led to more 



proposals and their experimental tests shows that Wigner's proof of linearity was not con- 
clusive; his assumption was questioned. I described where the proofs of linearity stood (24i] 
at the time of Weinberg's proposal. Soon thereafter, an idea was proposed that relativity 



requires quantum dynamics to be linear 
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22 , UM ■ This is discussed in Section IIVI 



II. LINEAR QUANTUM DYNAMICS 

The result of quantum dynamics is the time dependence of mean values for Hermitian 
operators representing physical quantities. This includes the time dependence of probabil- 
ities, which are mean values for projection operators. The result is the same whether it is 
obtained from the Schrodinger picture or the Heisenberg picture. The basic linearity to be 
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considered here is that the equations of motion for the mean values are hnear in that the 
time derivative of each mean value is a function of mean values that depends on the state 
in a linear way. This does not mean that the equations of motion for operators are linear. 
The time derivative of an operator may be a product of operators, so the time derivative 
of the mean value of the operator is the mean value of a product of operators. The basic 
linearity to be proved here is just that the time derivative of a mean value does not involve 
products of mean values. It depends on the state in a linear way. This is the case in ordinary 
linear quantum mechanics. Proposals for nonlinear quantum mechanics provide examples 
with equations of motion where time derivatives of mean values are functions of mean values 
... depend o„ .he .a.e . a nonlinea. wayfl Q . They involve p.odnc. of ,nean value. 
That happens because the Hamiltonian operator contains mean values; it depends on the 
state. 

The mean value (Q) for a Hermitian operator Q is Tr [Qp] where p is the density matrix 
that represents the state. In the Schrodinger picture, the time derivative of {Q) is 

dp 



(2.1) 



At each time, the equations of motion give time derivatives d{Q)/dt as functions of mean 
values at that time. The mean values at that time can be changed physically only by 
changing the state at that time. The basic linearity considered here is that each time 
derivative d{Q)/dt is a function of mean values that depends on the state in a linear way. 
The basic linearity fails if and only if there is a density matrix 

p = pp^ + {l-p)p2, (2.2) 

where pi and p2 are density matrices and p is a number between and 1, such that dp/dt 
is not the same as 



The basic linearity is that dp/dt and dp/dt are always the same. This means that density 
matrices are mapped linearly to density matrices in an interval of time. It is shown in 
Section Hm that with additional assumptions this leads to the linear form of the Schrodinger 
equation. 



Here is a proof that dp/dt and dp/dt must always be the same. It assumes the dynamics 
is for a system S that does not interact with anything else. Hence it assumes the dynamics 
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does not depend on anything that happens outside S, or on the state of any other system. 
In particular, this means the dynamics can depend on the state of S but not on any other 
property of the state of a larger system that contains S. It is also assumed that the dynamics 
must allow S to be described as part of a larger system. 

Suppose 5* is one of two separate systems 5* and R and that the state of the larger system 
of S and R combined is represented by the density matrix 

n = pp,\a){a\ + {l-p)p2\(3){(3\ (2.4) 

where |a) and \j3) are orthonormal vectors for R that do not depend on the time and, as 
before, pi and p2 are density matrices for S and p is a number between and 1. The reduced 
density matrix Tr^jII, which is the density matrix p for 5*, is described by Eq. ()2.2l ). The 
probability (P) for a proposition represented by a projection operator P for S is the sum of 
joint probabilities 

(P) = TTsR[P^=TTsB\P\a){a\^ + ^TsR[P\f3){m 

= {P\a){a\) + (P|/3)(/3|) = pTrs[Ppi] + (1 - p)Tr ^[Ppd. (2.5) 

Suppose a measurement is made on R that distinguishes the states represented by \a) and 
1/3). The probability is p that the result is \a) and 1 — p that the result is If the result 
is la), the probability for the proposition represented by P is Tr5[Ppi], and if the result is 
1/3), the probability for P is Tr5[Pp2]- This can be verified experimentally by repeating the 
process of preparing the state represented by 11, measuring to distinguish the states of R, 
and testing various propositions for S. The density matrices pi and p2 describe physically 
distinct possibilities. The times of the events can be changed as long as the measurement on 
R is early enough. The time dependence of pi and p2 will account for changes in the results. 
The time derivative of the probability for the proposition represented by P is TTs[Pdpi/dt\ 
or TTs[Pdp2/dt] depending on the result of the measurement on R. Altogether, with the 
probabilities for both results being considered, the time derivative of the probability for P 
is TTs[Pdp/dt\. 

The probability (P) for the proposition represented by P is also 

(P) = Tr5,j[Pn] = Tr5[PTr^n] = Tr5[Pp]. (2.6) 

Its time derivative is TTs[Pdp/dt\. This is always correct. It may be the only possibility at 
hand. For example, suppose the state of S and R combined is represented by the density 
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matrix 

n = pp\a){a\ + (1 - p)p\(3){f3\ = p[p\a){a\ + (1 - (2.7) 

The dynamics for S must be the same for 11 as for 11. 

The time derivative of the probabihty for the proposition represented by P is always 
TTs[Pdp/dt]. There are situations where it also must be Tr s[Pdp/dt] to fit observations of 
events in a larger system. The dynamics in 5* can not depend on the situation of S* in a larger 
system. Therefore TTs[Pdp/dt] and TTs[Pdp/dt] must be the same. From that equality for 
various projection operators P, we conclude that dp/dt and dp/dt must be the same. 

Nothing in this proof requires quantum mechanics. Quantum entanglement is not used. 
The correlations involved are classical. The substance would be the same in classical me- 
chanics; only the language would be different. 



III. THE SCHRODINGER EQUATION 

The result of Section |n] brings us within sight of a proof that the dynamics can be 
described by changes of state vectors in time determined by a linear Schrodinger equation. 
The classic proof assumes that in an interval of time the dynamics maps the set of all pure 
states one-to-one onto itself. The first step in the proof is Wigner's theorem that if the 
absolute values of inner products of state vectors do not change, the change of state vectors 
can be made with an operator that is either linear or antilinearfl 0. The product of two 
antilinear operators is linear, so if the change of state vectors from one time to another can 



be made in two steps, it must be linear j26|. If the change of states does not depend on the 
time when the change begins, and if the change of probabilities in time is continuous, then 
the state vectors as functions of time can be obtained with a continuous one-parameter group 
of unitary operatorsjs^ 22, 0|. They satisfy a linear Schrodinger equation; the Hermitian 



generator is the linear Hamiltonian operator [ 

If the dynamics applies to all states, it must map every density matrix p to a density 
matrix p' in an interval of time. If density matrices pi and p2 are mapped to p[ and p'2, 
then what was proved in Section |n] implies that the density matrix described by Eq. ()2.2I ) 
is mapped to 

p' = pp[ + (l-p)pl (3.1) 
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The map is linear for density matrices. It has been known for some time that this can be 
used to prove the assumptions of Winner's theorem and open the way to the hnear form of 
the Schrodinger equation j^, 13- 

Linear maps of density matrices can also be used to describe processes where pure states 
are mapped to mixed states, different states are mapped to the same state, the map is not 
onto all states or, generally, the map has no inverse that applies to all states. An assumption 
is needed to separate these processes from dynamics described by the Schrodinger equation, 
which has an inverse for all states. 

Suppose that in an interval of time the dynamics maps the set of all density matrices 
one-to-one onto itself. The result of Section |n] is assumed, that the map is linear for density 
matrices, as described by Eq. (jH.lj) for a density matrix described by Eq. ()2.2|l . Then the map 
has an inverse and the inverse map is linear; the proof that the inverse of a linear operator 



is linear[26|. Theorem 7.1] applies with attention restricted to density matrices. Pure states 
are mapped to pure states: if (J^ and p'2 in Eq. ()3.Hl are distinct, so are pi and p2 in Eq. ()2.2|l : 
thus if p' is for a mixed state, so is p. The inverse map also takes pure states to pure states, 
so the set of all pure states is mapped one-to-one onto itself. 

For each vector of length 1, let be a vector of length 1 such that dV') ("^D' is 
I'?/'') ('?/'' I . For each density matrix p there are orthonormal vectors and positive numbers 

whose sum Yl,jVi is 1 such that 

P = Ep.IV^.)(^.I- (3-2) 



The linearity implies that 



Y.vM'm\- (3-3) 



Since ^\^'M)'^ is 1, 



Tr[(p')'] = Y.vMm^'^)\^ > Efe)' = Tr[/]. (3.4) 

jk j 



The same result for the inverse implies that 

Tr[(p'f] = Tr[p2]. (3.5) 

Let and |0) be vectors of length 1 and let 

P = ^IV^)(^I + ^I0)(0I- (3.6) 
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Then 



Tr[p2] 



1 



1 



1(^10)1 



2 



(3.7) 



2 



+ 



2 



and 



p' = \w)m + \w){^' 



(3.8) 



so Eq. ()3.5I ) implies that 



m')? = \m)?- 



(3.9) 



Absolute values of inner products of state vectors are not changed. The assumptions of 
Wigner's theorem have been proved. The derivation of the linear form of the Schrodinger 
equation can proceed along the classic route. 

An alternative is to assume that it is the set of all pure states that is mapped one-to-one 
onto itself, instead of the set of all density matrices. The result of Section |n] is assumed as 
before, so the dynamics defines a map of density matrices to density matrices that is linear 
for density matrices. Then the inequality ()3.4|) is obtained as before. If it is assumed that 
the entropy does not decrease|3], then 



This implies Eqs. fl3.5p and ()3.9|) . An alternative to assuming that entropy does not decrease 
is to assume that the probabilities in mixtures of orthogonal pure states do not change [2^ . 



Questions about previous proofs helped motivate this work. The idea that relativity 
might require quantum dynamics to be linear was clearly expressed in an example j3] where 
nonlinear dynamics makes a dp/dt defined by Eq. ()2.3|l different from dp/dt and lets it be 
changed from outside the system. The system S is in an entangled state with a separate 
system R so that one of two possible results of a particular measurement in R implies 
that probabilities of propositions for S are described by the density matrix pi, and the other 
possible result from R implies that probabilities for 5* are described by the density matrix p2- 
Let p and 1 —p be the probabilities of the two possible results from R. The dp/dt of Eq. ()2.3|l 
is not the same as dp/dt. A different measurement in R yields different density matrices pi 
and P2 and a different dp/dt. If dp/dt describes the time dependence of probabilities for 5*, 



Tr[(p')^]<Tr[p2]. 



(3.10) 



IV. DISCUSSION 



8 



then observations in 5* can determine which measurement was made in R. A signal can be 
sent from R to 5*. Since there is no restriction on the locations of R and S, the signal can 



be faster than light. This assumes that dp/dt, not dp/dt, describes the time dependence of 
probabilities for S. This assumption can be questioned |3, 

If a signal can not travel from i? to 5* at the speed of light in the time between the 
measurement in R and observations in S, then moving observers can disagree about whether 
the observations in S happen after or before the measurement in R. The probabilities for 
the results of the measurement in R are calculated from the state of the larger system of 
S and R combined, without consideration of what happens in S. Probabilities for S can 
be calculated the same way, from the density matrix p for 5* that is the reduced density 
matrix obtained by taking the trace for R of the density matrix for S and R combined, and 
is related to pi, p2, p and 1 — p by Eq. ()2.2|) . The state of 5* and the state of R both are 



prepared when the state of S and R combined is prepared. Can we not just say 



19, 



33, 



that the time dependence of probabilities for S is described by dp/dt, not dp/dt, in this 



situation? Then there would be no signal faster than light. Use of dp/dt could be reserved 
for mixtures where dpi/dt and dp2/dt have to be calculated separately because the states 
represented by pi and p2 are definitely prepared separately^^. Should we not be prepared 
to refine our rules of interpretation of quantum mechanics just this way if experiments show 
evidence of nonlinear quantum dynamics? 

If there is time for a signal to travel from i? to S" at the speed of light and report the result 
of the measurement in R before the observations are made in S, we can say that a state of 
S represented by pi or p2 is definitely prepared before the observations, so that altogether, 
when the probabilities p and 1—p for both possibilities are considered, the time dependence 



of probabilities for 5* should be described by dp/dt. All of this would be simpler in the 
context of nonrelativistic mechanics where moving observers agree on the order of events in 
time and there is no limit on the speed of signals. Could relativity actually provide a niche 
where nonlinear quantum dynamics can exist? If there were nonlinear quantum dynamics. 



there would have to be an abrupt change from dp/dt to dp/dt when gradual changes in 
distances and times reach the point where a signal can no longer travel from i? to S* at 
the speed of light in the time between the measurement and observations [ssl . This makes 
nonlinear quantum dynamics look unreasonable, largely unattractive for investment of time. 
Does it prove that nonlinear quantum dynamics is impossible? That it does not can be 
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demonstrated with a fanciful construction. j3f 

This example uses Weinberg's nonlinear quantum mechanics. When the idea that relativ- 



ity implies linearity was presented more generally 
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the proofs assumed that if a map 



of density matrices to density matrices is defined by the dynamics, it must be linear in mix 



tures of pure states. This does not hold for Weinberg's nonlinear quantum mechanics 
formulated in terms of density matrices jiol. l2^ . so these proofs do not effectively address 
Weinberg's proposal; it is simply ruled out by that assumption at the start. The same 
assumption was made to present the idea that nondecreasing entropy requires hnearityjl^. 
Using entropy together with basic linearity as in Section IIIII simplifies the proof using en- 
tropy and makes it clear that the only additional assumption that is needed j2J] is that pure 
states are mapped one-to-one onto pure states. 



stops when a linear 
25l | goes on to a proof 



The first presentation of the idea that relativity requires linearity 
map of density matrices is obtained. The more recent presentation 
that the map of density matrices is completely positive, so that it "can be realized ... by 
a linear and unitary evolution on a larger Hilbert space." This would include dynamics 
that is not described by a Schrodinger equation but is for a subsystem of a larger system 
where the dynamics is described by a Schrodinger equation. It would bring in cases where 
the considered system interacts with other parts of the larger system. We know now that 
evolution in a subsystem caused by linear unitary evolution in a larger system is described 
by linear maps of density matrices that generally are not completely positive and act in 
limited domains [2^. Complete positivity generally will not provide a link to the source of 
the dynamics. Including subsystem dynamics will be complicated by the need to consider 
domains. It is simpler to consider only the given system and make the assumptions needed 
to link up with Wigner's theorem as in Section IIIII 
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